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O . We review the theory of the determination of the CKM elements | Vcb \ and | Vub \ ■ Particular 



attention is paid to the determination of \Vcb\ through inclusive semileptonic B decays to 
charm using a moment analysis, since this has shown most progress recently. A precise 
. method for the determination of \Vub\ via exclusive decays is described. 



> 

2h : 1 Introduction 

^ ■ The Wolfenstein parametrization of the CKM matrix, 

: V^cKM = Vcs Vcb\^\ -Xil+iA^X^r) 1 - +0{\'') (1) 

Ph: \Vtd Vs VbJ \AX\l-p-zr]) -A\\l + iX'r]) 1 / 

D ! shows explicitly that it is determined by four parameters. Ki3 decays determine A to 1% 
[ accuracy!!]. The other three parameters are the subject of this talk. As we will see, \Vcb\ = AX'^ 
is now determined to nearly 1% accuracy from B — > Xcii^e- The magnitude of the remaining 
^ ■ complex parameter, p — ifj = Vub/iVusVcb), is fixed from B Xuii^e, but is known far less 
accurately. 

Joining the points 0, 1 and J) + ifj in the complex plane yields a unitarity triangle. The sides 
have length 1, \/^^TW = iKbl/lK^Kbl and - p)^ + rj^ = \Vtd\/\VusVcb\. To determine 
the triangle one needs not just the magnitude but also the phase of p — ir]. Alternatively, the 
length of the third side determines the triangle, up to a two-fold ambiguity. This third side can 
be obtained from b ^ d'y decays (eg, B p'y) or from the B^ — B^ mixing parameter AM^. 
In the standard model, both of these processes are dominated by a virtual top quark exchange 
and depend directly on V^^V^l- The theory of these processes belongs in a different session in 
this conference 12 El- 

The bulk of the talk is devoted to the determination of \ Vcb\ from inclusive decays B — > Xcii^e, 
since here is where we have seen tremendous progress in the last year. Next we briefly revisit 
|Vcfe| from exclusive decays, B — »• D*ii'£ and B —>■ Diue. We then turn to \Vub\ and mirroring 
the discussion above we first present its determination from inclusive decays, B —>■ Xuii^e- This 
is followed by a review of the extraction of \Vub\ from exclusives B —>■ 7r{p)ii'£. The limitted 
precission in the determination of \Vub\ afforded by the standard methods calls for innovation 
in this area. We describe a new method that uses a combination of radiative and semileptonic 
B and D decays. 
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2 \Vcb\ from Inclusive decays B Xcii^i 

In QCD the differential width for B — > XcCvi is a function of \Vcb\i the strong couphng as and 
the quark masses m{, and mc. Only the dependence on \ Vcb\ is trivial: 

dT{B ^ X^y,) 

-^j^^ = \Vcb\ f{x,y;mb,c,as), (2) 

where x, y are kinematic variables. Even if we could compute the function / from first principles, 
we would still need as input precise values of the three parameters, nib, f^c and as- While as 
is accurately known J^], the quark masses are not. For a determination of Vcb one can use the 
measurement of the doubly differential decay rate to fit also the value of the masses. 

Now, Vcb enters only the overall normalization, while rrib and rric determine the "shape" 
[i.e., the functional form) of the decay distribution. So one can isolate the determination of the 
masses by fitting to shape parameters, that is, various moments of the distribution. In all cases 
the moments are defined with a charged lepton energy cut, an experimental necessity turned 
into a theoretical tool. From the charged lepton integrals with an energy cut 

Rn{E,^u M)= ! {Ee - M)"-^ dEe , (3) 



define moments 



,pn\ _ Rn{Ecut,0) ..X 

\Ee /iJcut = s-TB FTV • \V 



RoiEcut, 0) 
Central moments are also defined: 

/(TP RnjEcut, {Ee)Ecut) /r\ 

\\Ee - {Et)E,^J )e,^, = ^-TT^ ■ [o) 

-Kol^-c/cut, Uj 

Also useful are moments of the hadron invariant mass, = {pb — Pe — PvY — (Pb — qY-, 
defined by 

)E,ut = r dr 1 2 ■ 



Central moments ((m|^ — m|))"), about = [{niD + 3mD.)/4]^ are also frequently considered. 
Somewhat orthogonal information can be garnered from moments of the B Xg'y spectrum, 
which depend on m& and as, and only very weakly on rric- They are defined for the spectrum 
of the photon energy in the B rest frame, E-y-. 

I rpn\ ^Ectii-^'^^ dEj *^-^7 , , 

= f ji^dE 

Ecut dE^ 7 

The variance, {E^) — (E^)"^, is often used, but higher moments are not, because they are very 
sensitive to the boost of the B meson and to details of shape function. 
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2.1 Heavy Quark/Mass Expansion 

At present we don't know how to perform the calculation described above. Instead a systematic 
approximation to QCD is made by expanding in inverse powers of the heavy mass, rrib, and 
performing an operator product expansion[S]. The price we pay is that new dimensionful 
parameters are introduced. At n-th order in the expansion the new parameters are roughly 
of size (tIqcd)'^, giving rise to an expansion in powers of {Aqcb / n^'b) ■ These new parameters 
account for our ignorance of non-perturbative dynamics in QCD. 
There are two approaches to the expansion in the literature: 

• The HQET/OPE, in which amplitudes are given in terms of mass independent heavy 
meson states. If one chooses to treat the iI)''*^-meson as heavy then this approach involves 
an expansion in l/rric (in addition to l/rrib). The advantage is that rric can be fixed in 
terms of ruh using the accurately determined — Mb- 

• The HME, or Heavy Mass Expansion, in which heavy meson states are not expanded 
(in neither rric nor mf,). No l/rric expansion is necessary. However, rric is a parameter 
determined from the fit. 

Three groups have performed a moment analysis. The groups of M. Battaglia et and 
of Gambino and Uraltsev[Zj use the HME, while C. Bauer et a/[Hj consider and compare both 
expansions. In all cases the expansion is up to (and including) terms of 0{Aqcj))^, and known 
perturbative corrections are included. More on this below. 

2.2 Parameter Counting 

It may come as a surprise that the number of parameters introduced by the two approaches is 
the same. The naive guess that the HQET/OPE with a l/rric expansion ought to introduce 
more parameters because it expands in l/rric too is incorrect. The additional expansion gives 
new information, for example, that to leading order the B and D meson states form heavy 
flavor symmetry doublets. 

Altogether, there are six parameters (in addition to as) that enter the moment analysis 
in either approach. With a l/rric expansion rric is fixed using Mb — M^, and the parameters 
A2 ~ (tIqcd)^ and p2 ~ (^qcd)^ are fixed from Mb* — Mb and M^)* — Mb- One is left with m^, 
Ai ~ (/Iqcd)^, and four parameters of order (/Iqcd)^, namely pi and three linear combinations 
of four time ordered products, 71-4. Now one has the issue of accuracy of expanding in l/ntc- 
The OPE is an expansion in 1/mfe. It is only the computation of ttIc in terms of meson masses 
that requires a l/rric expansion. Since this is done to 0{AQQYi/^c)'^ the fractional error in 
rric is 0{AQQj)/m()'^- But mc first enters the rate at 0{mc/mhY , so the error introduced is 

C?(4cdM>D- 

If the states are not expanded in 1/m (so no l/rUc expansion is performed) then no time 
ordered products appear. However, rric, ^2 and p2 cannot be solved in terms of physical meson 
masses, and they come in as parameters. Comparing with the counting above, the three 
combinations of time order products are replaced here by the three parameters rric, ^2 and p2- 
Strictly speaking the parameters should be distinguished from the ones above and denoted by 
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different symbols since in the HME the non-perturbative parameters have imphcit dependence 
on the heavy quark mass. Bauer et al suggest that the fit to moments may be better behaved 
when the 1/mc expansion is performed, since more of the parameters are of higher order in 

^QCD- 

In either approach one needs to properly define quark masses. A good definition of the quark 
mass will give better convergence of the perturbatively calculable (in a^) Wilson coefficients of 
the OPE. Many definitions are available^. Threshold mass definitions, like the IS and PS, give 
better convergence than the pole, MS or kinetic mass schemes. Battaglia et al and Gambino 
and Uraltsev use the kinetic mass scheme, while Bauer et al use all and compare. This latter 
group points out a problem with the kinetic mass. Since this scheme requires a rather low hard 
cut-off /i, when there are two heavy quark masses, mi, and nic one introduces two distinct hard 
cut-offs, Hb and /ic. While it is customary to use /ib = /ic, there is nothing in principle which 
requires one to use equal cut-offs for the two masses. Bauer et al consider the effect of changing 
yUfe^c so that /ift 7^ /ic, and find rather large variations in the determination of \Vcb\- The effect 
is due, presumably, to the next-order corrections in ^^(/ib.c), which are rather large given that 
libfi must be taken to be of the order of the fairly low scale ~ 1 GeV. 

2.3 Fit results 

The BaBar Collaboration has performed a fit to the moments calculated by Gambino and 
Uraltsev. The electron energy moments are computed to order /3oa^, cXs^AqqdY ^-^d (/Iqcd)^, 
while the hadronic mass moments are computed to order a^, ^^(^qcd)^ and (tIqcd)^- To 
account for the missing order PqoI effects, the analysis uses a different value of as for the 
hadronic mass moments than for the lepton energy moments. This educated guess is somewhat 
ad hoc and introduces an unknown uncertainty. The result of the fit is presented in this 
conference elsewhere The fit includes non- integral moments of m^, which are theoretically 
less well behaved and have larger uncertainties. 

The group of Bauer et al performs a global fit using data on semileptonic and radiative 
decays from the BABAR, BELLE, CDF, CLEO, DELPHI collaborations. They compare the 
two approaches, expanding or not in m-c, and the different mass schemes. Fig. Q shows the 
convergence in of the computation in different schemes. The dotted (red) dashed (blue) and 
solid (black) ellipses denote the result of the fit when retaining terms of order a^, and /5oa^, 
respectively. 

Figure 121 shows the result of the fit, for several moments as a function of E^nt- The data points 
are superimposed: square (blue), triangles (green) and filled circles (black) are from BABAR, 
CLEO and BELLE, respectively. The half integer hadronic moments are shown, but not used 
in the fit for |V"cf,|. The shaded region indicates the theory error, which has been estimated 
from the size of the highest order retained in the double expansion in as and Aq^cTi/^b- The fit 
shown here is for the HQET/OPE approach (with 1/mc expansion) in the IS mass scheme. The 
graphs labeled BR and Ml show the branching fraction and the first lepton energy moment, 
while M2 and M3 show the second and third lepton energy central moments. The remarkable 
agreement between theory and experiment indicates that the expansion is working very well 
and, in particular, that local duality works extremely well for these decays. The result of the 
fit for \Vcb\ and irib is shown in Fig. El It is quite remarkable that the fitted value of irib is 
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Figure 1: Convergence in of the fit, in different schemes. The dotted (red) dashed (blue) 
and sohd (black) ellipses denote the result of the fit when retaining terms of order a^, as and 
Poal, respectively. 




Figure 2: Fit of moments as a function of -E'cut- The data points are from BABAR (squares), 
CLEO (triangles) and BELLE (circles). 
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Figure 3: Result of fit for IV^^I and m^. The lower graphs do not include theory errors and can 
be used to see the present theoretical limit of the extraction of parameters. The ellipses denote 
regions of Ax^ = 1 (black and yellow), A^^ = 4 (blue and green). For the ellipses in black 
and blue the parameters of order {AqcbY cire restricted to vary within ±(500 MeV)^ while no 
such restriction is used in the fit shown in yellow/green. The red cross shows the current PDG 
value of \Vcb\ and the IS rrib determined from a fit to bottomonium(lIj. 

comparable to that obtained from masses of bottomonium and with comparable precision JT]. 

3 \Vcb\ from Exclusives: B D*iu£ and B Diu£ 

The rates are given by the HQET-inspired parametrizations 



dw 



dT{B DtV) 
dw 



487r^ 



1 + 



Aw 1 — 2wr^: + rz 



l + w 



^i^r^{l + rf{w'-lf''\V^b?r\w) 



in terms of T and which are simply combination of form factors oiV — A charged currents. 



Here w = pb'Pd /MbM}^\ and r^*) = M}^' /Mb- At lowest order in HQET J^(l) = = 1, 

and Luke's Theorem insures the corrections are small, — 1 = 0{AQCD/mc)^. 

In order to determine \Vcb\ we need an accurate theory of a measurement of the 

rate as a function of w and an extrapolation of this measurement to w = 1. Analyticity and 
unitarity tightly constrain the functional form of ^*(w)[T^. Taylor expanding about w = 1, 



J^*H = + p^iw - 1) + c{w - 1)2 + . . .] 



(8) 



the unitarity constraint gives a relation between the slope, p^, and curvature, c^^. Sum rules 
give lower bounds on these parameters ^3], but the bounds are not experimentally significant. 
It is interesting to note that at leading order in HQET the slope and curvature parameters are 
equal for J-'{w) and J-'^{w). Corrections have been estimated, — ~ 0.19, to be compared 
with the experimental fit, — p^^ ^ —0.22 ± 0.20^2I- A more precise determination of the 
slopes and curvatures would help (and may tax) the theoretical understanding of sub-leading 
form factors in HQET. 

The main limitation in determining \Vcb\ comes from theoretical uncertainties in J?%=(1) — 1. 
The situation has not changed since the BaBar book was pubhshed, JF^(l) = 0.91 ± 0.04[Tnj. 
The current lattice determination has comparable errors, ^*(1) = 0.913 1q;o^7 ^aoso CS]- For a 
comparison of the determinations of \Vcb\ from inclusive and exclusive decays see Ref. jTUI in 
these proceedings. 

4 \Vub\ from inclusive B Xuii'i 

The theory of inclusive b — >■ uiu processes is the same as in 6 — ciu. However, a straightforward 
application of theory is not possible at the moment since experimentally the full spectrum is 
not available. The problem is that the signal is swamped hj B ^ Xcii^e- The experimental 
solution: impose kinematic cuts to suppress the b —>■ c background. By imposing either, 
mx < rriD, > {niB — m£,Y or > {m\ — m|,)/2mB one restricts the measurement to events 
which are not kinematically accessible to 6 — > c decays. However, the experimental solution is 
largely incompatible with theory [20] ■ The endpoint spectra m B ^ X^lve, is given in terms of 
a non-perturbative "shape function," f{x). One approach is to measure f{x) in B and 
use it to eliminate uncertainties in the determination of \Vub\ from the B Xud-Vi analysis [T??j. 
However, the universality of the shape function is violated at order l/ruh, and the size of these 
corrections is uncertain 

The sensitivity to the shape function is least with the cut on the lepton invariant mass, 
> {niB — moYi but the cut picks up a small region of the Dalitz plot. The hadronic mass 
cut includes a large fraction of the Dalitz plot, but is more sensitive to f{x). For now, the best 
results are obtained by combining both cuts, trying to maximize the rate while keeping depen- 
dence on f{x) to a minimum. Theoretical estimates of the uncertainty in the determination of 
\Vub\ from this method are in the 10% range [22] • 

5 \Vub\ from Exclusive Decays 

The determination of \Vub\ from exclusive decays requires a priori knowledge of the form factors 
for i? ^ TT and/or B ^ p, e.g., 

{7r{p^)\h,q\BiPB)) = (^(PB+P.), - "^^^^q^j F+{q') + ~ g.i^o(g^) • (9) 

HQET does not fix normalization at zero recoil as it does in i? — > D^*^^. Analyticity and 
unitarity do constrain the functional form, which helps, e.g., to interpolate lattice results. 
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Figure 4: Form factors for B Tciu computed by simulations of lattice QCDfTBj. 



An interpolation of lattice results is shown in Fig. |3J One can use this to determine \Vub\, 
but the errors are large. Belle finds, from B iriu restricted to > 16GeV^, \Vub\ = 
(3.87 ±0.70 ±0.221°:^^) x lO'^ and |Kb| = (4.73 ± 0.85 ± 0.27l[5;^^) x 10-^ using lattice results 
from FNAL'04 and HPQCD, respectively[2ni- The theory errors, ~ ±20%, are not expected 
to be significantly reduced soon. There has also been some effort to understand the precision 
with which \Vtd\ can be determined from B ^ p'y decays pij. 

New ideas are needed to reduce the error on \Vub\ from exclusive decays to the sub-10% 
level, hopefuly to a few percent. One recent proposal is to use double ratios [2^1 to eliminate 
hadronic uncertainties. A double-ratio is a ratio of ratios of quantities that are fixed by two 
symmetries. For example, 

R^ = ^ = ^j^ = l + o(mJ^-l-)) (10) 
JdJId Jb/Jd V mbJJ 

The ratios /bs/ fs and fnJfD are unity in the fiavor-S'f/(3) limit, while the ratios /bs/ fn^ 
and Jb / fo are fixed to a/ rric/mb by heavy-quark symmetry. Hence Ri — 1 vanishes either as 
ms — * or as mq — > oo, as indicated. 

To use a double-ratio in a precise determination of |\4fe/\46Vts|, one measures, for 4m^ < 
(f ^ 9max = {rnB-mp^K*f 

dr(5^irW-)/dg2 |\4,V;*|2 a2 iVeff(g2) ^Jiff-^' (g2)|2- ^ ^ 

Here Hx are helicity amplitudes, and A'efr(g^) is a calculable function that incorporates the 
effects of penguin diagrams. In addition one must measure decay spectra for D piu and 
D — > K*ii> and express all rates as functions of y = Ey/mv {V = p,K*). Then one uses 
double-ratio magic. Let 

. EJgf^%)P . E,\Hr''iy)\' 
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then 

RB^v{y) = RdMv) (l + Oim.i- -))) (13) 

V rric rrib / 

can be used to eliminate the unknown hehcity amphtudes from the analysis. 

For the calculation of the factor iVefr(g^) = \Cf + ^Cfl"^ + \Cio\^ + 0{A/mb), one may 
include the effects of penguin diagrams through an OPE. This is indicated above through the 
superscript "eff" (Cy^g C'rg)- In applying the OPE, there is no need to assume that the 
four quark operators factorize into the product of two currents. However, the OPE is used in 
the timelike region, so this assumes local duality. The situation is under fairly good control, 
however, since the penguin contribution is very small compared to the leading contribution from 
the contact terms. Although the contribution is small, it is crucial to include it in order to 
minimize the scale dependence of N^siq'^)- The uncertainty in the scale dependence in Neg{q'^), 
which should be an RG-invariant, is only a couple of percent in a NNLO calculation. 

Although a comprehensive study of the corrections has not been performed, preliminary 
estimates indicate that the theoretcial error in the determination of \Vub\ could well be bellow 
10% with this method. For example, the deviation of Ri from 1 is about 2% in model and 
lattice calculations |27j. Similar results are obtained for the double ratio of form factors for 
B/D K*/p^. 
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